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Quantum fluctuations created during inflation can account for the observed matter distribution 
in the linear regime if the universe has two components of dark matter, one which is cold and 
collisionless, and the other which is hot and free streams on small scales. However, this free streaming 
property of the hot component prevents early structure formation, and since objects, such as damped 
Lyman-a systems, have been observed at high redshift, it is necessary to produce more power on 
small scales. Here, we show that the situation can be improved substantially in models where cosmic 
strings are formed at the end of inflation, and in which both inflation and strings participate in the 
generation of structure. 
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During the next decade various high precision mea- 
surements will enhance our understanding of the universe 
immeasurably. Two satellites (MAP and PLANCK) will 
map the cosmic microwave background (CMB) over a 
wide range of scales and ground-based redshift surveys 
(for example, Sloan Digital Sky Survey and 2Df) should 
measure the clustering of luminous matter. Using these 
complementary measurements — as well as many others 
(for example, those of Type la supernovae and light el- 
ement abundances) — it is hoped that robust estimates 
can be made for a large number of cosmological parame- 
ters. However, underlying this ambitious enterprise is the 
origin of the primordial fluctuations, which is usually as- 
sumed to be a near scale invariant spectrum of adiabatic 
density perturbations created by quantum effects during 
some de-Sitter phase associated with cosmic inflation. 

Under various assumptions as to the matter content of 
the universe three broad classes of models based on infla- 
tion have been suggested to fit the current observations 
of galaxy clustering in the linear regime (A > 30ft.~^Mpc) 
when normalized to the amplitude of anisotropics in 
the CMB detected by the COsmic Background Explorer 
(COBE) satellite (see, for example, ref. Q). Each of 
these models is based primarily on a universe whose mat- 
ter density is dominated by cold dark matter (CDM) par- 
ticles {flm ^ ^c), such as neutralinos or axions, with a 
much smaller component of baryons as predicted by Big 
Bang Nucleosynthesis fib ~ 0.05 — 0.1. At present, pos- 
sibly the most popular of these classes of models have a 
matter density which is less than critical (flm ~ 0.3—0.5), 
with in one class (ACDM) the shortfall from critical being 
made up by a non-zero cosmological constant (r^A+^m — 
1) and the other (OCDM) having an open topology. The 
final class of models (CHDM) M have fim = 1 with 
an extra relativistic or hot dark matter (HDM) compo- 
nent {Qi, w 0.2 — 0.3), such as massive neutrinos, which, 
since they free-stream, reduce the amount of power on 



intermediate and small scales, allowing the model to fit 
the shape of the observed power spectrum in the linear 
regime. Reconciling the predictions of this cold plus hot 
dark matter model with observations at high redshift is 
the subject of this letter. 

The relevant observations concerning us are those of 
damped Lyman-a (DLYA) systems Q which are thought 
to probe structure as it was for z > 2 in the linear regime, 
and on scales smaller than otherwise possible. There 
seems to be a conflict between the measured neutral gas 
fraction, particularly at very high redshifts (z « 4), and 
the predictions made for CHDM models MM, which — 
due to the free-streaming property of the HDM compo- 
nent — tend to lack of power on such small scales. This 
is one of the main reasons why, in spite of the superior 
achievements of CHDM models on linear scales, OCDM 
and ACDM models are now preferred. 

For many years the most popular alternative to the 
standard adiabatic scenarios discussed above has been to 
form structure by accretion of matter onto topological 
defects (for a review, see ref. Q). After much recent 
work [|7|-[10[, it appears that, at least when Om = 1, these 
models appear to suffer from the converse problem; they 
have a tendency to produce insufficient power on large 
scales to match the current observations when normalized 
to COBE, but are known to produce copious amounts of 
small scale power [p|-pX)|, even in scenarios with a HDM 
component uM. This is particularly true in the case of 
cosmic strings and one may wonder if the complementary 
failings of CHDM based on inflation and string models 
could be erased in a mixed string plus inflation CHDM 
model. 

The possibility that cosmic structure could be due to 
both strings and inflation |12| is a very real one which 
has been the topic of much recent research 13-lq|, since 



it could arise in a number of well motivated inflation- 
ary scenarios pq-EOl . The details of this paper are only 



weakly sensitive to the specifics of the inflationary model 
and then only in an essentially predictable way. However, 
to give an example, let us consider D-term inflation |1^] 
(see ref. |2^ for a review). These models make use of 
supersymmetry in order to provide cosmology with a flat 
potential, as required by inflation, without the need for 
a finely tuned coupling constant. During inflation, the 
inflaton moves along this flat direction, which is a direct 
result of the U(l) symmetries in the model, and at the 
end of inflation the underlying symmetry becomes bro- 
ken, producing cosmic strings via the Kibble mechanism 
during the subsequent phase transition. 

In any string plus inflation scenario the large angle 
CMB anisotropics, which are normalized to COBE, are 
due both to strings and inflation with a model dependent 
ratio between the two; it being a non-trivial function of 
the number of e-foldings during inflation, the inflationary 
energy scale and that for the symmetry breaking phase 
transition which produced the strings. Given a specific 
model it would be possible to relate this to parameters of 
the model as was done in refs. ||l^,|l^ . In both cases it was 
shown that any value of this ratio was possible for sen- 
sible model parameters and therefore for the purposes of 
this paper we will leave it arbitrary. In particular we will 
use the notation of ref. |Tq] , where the power spectra were 
added as, for example, P{k) = aP'"'^{k) + {1- a)P"*''(fc), 
where < a < 1 and P'"^(fc), P^^''{k) are individually 
computed spectra for the inflation and string models nor- 
malized to COBE. The specific spectra that we shall use 
have been computed using CMBFAST Q and in the 
string case the model for the defect stress energy used 
in refs. PJl5||. This has been shown to represent many 
of the features of a realistic string network, although it 
may under estimate the amount of small-scale power pro- 
duced 1^,0. Any extra small-scale power will further 
improve the situation relative the observations being dis- 
cussed here. 

The neutral gas fraction in DLYA systems has been 
computed using a variety of techniques for the stan- 
dard CHDM scenario based on inflation. Initially [Q|^] 
the Press-Schechter (PS) approximatio n p3^ w as used to 
make an estimate and more recently |24|-|27f| hydrody- 
namical simulations have improved upon this, although 
the basic picture has remained the same. In attempt- 
ing to estimate the size of the effect of including a string 
component one is faced with a number of obstacles to 
making a quantitative prediction, mainly related to the 
non-Gaussian nature of the perturbations involved. This 
makes creating a realization of initial conditions for a 
hydrodynamic simulation almost impossible without an 
accurate simulation of the defect network and of course 
one of the basic premises of the PS formalism is that the 
fluctuations initially form a Gaussian random field. It has 
been suggested 28 1 that a simple generalization of the PS 
approximation can be made by replacing the Gaussian 
probability distribution function (PDF) with that com- 



puted from a simulation and, indeed, it has been shown 
that such a modified theory yields good predictions for 
particular PDFs using TV-body simulations ||2^. A PDF 
has been computed for the pure string component pQ 30 
showing that it has positive skewness. In the mixed sce- 
narios considered here, the total PDF is the convolution 
of the individual components, and so the total cumulants 
are the sum of the individual cumulants. Hence, the com- 
bined PDF will also have positive skewness. 

Here, we will repeat the standard PS calculation as- 
suming Gaussian initial conditions following closely the 
calculation of ref. [Q . Since the measured skewness for 
the cosmic strings PDF is positive, one might expect that 
this calculation would act as lower bound on the gas frac- 
tion created in a given scenario. However, as we shall see 
the effects of non-Gaussianity are subtle and this is not 
always the case. To illustrate this we will use a simple 
non-Gaussian distribution with a PDF which is a log- 
normal distribution, as suggested in ref. ||29| . 

DLYA systems are observed as wide absorption troughs 
in the spectra of distant quasars. These troughs in- 
dicate that the line of sight to the quasar intersects a 
region of neutral hydrogen (HI) with a column density 
> lO'^'^cm^'^. For high redshift objects {z > 2) there are 
hints, that these absorption lines are produced in turbu- 



lent protospheroids |32|, which are the natural progen- 
itors of galaxies. The measured abundance of damped 



Lyman-a systems at redshift z = A |33|, implies that the 
neutral gas density f2gas in units of the critical density is 
^gas (2 ~ 4) ft, = (9.3 ± 3.8) X 10^*, in a universe which 
has Qrn — 1, where the Hubble constant is parametrized 
in the usual way Hq = 100/ikmsec^^ Mpc^^. 

The neutral gas fraction in DLYA systems is given by 
^gas = /ffl^berfc ((5c/v2cr_R) , which is the result of the 
integration over the tail of the underlying Gaussian PDF. 
In this expression < /hi < 1 is the fraction of hydro- 
gen which is neutral, flh is the average baryon density 
relative to critical, <5c(~ 1.69) is the density contrast re- 
quired for collapse to take place and a^ is the rms density 
fluctuation in a ball with size of the DLYA systems. To 
estimate this scale we assume that the systems are pro- 
tospheriods which are about to collapse to rotationally 
supported gaseous discs. In this spherical collapse sce- 
nario iQ the comoving size of the system is related to 
the circular velocity, T4, by 



R = 86/i"^kpc 



V, 



50km sec" 



1/2 



(1) 



A typical velocity is Vc ~ 50kmsec~^, which corresponds 
to the minimal mass {^ IO^^/i^^A/q) needed for spherical 
collapse to take place p5[. The rms fluctuations on this 
scale, an, are computed in the standard way by integra- 
tion over a window function corresponding to a spherical 
top-hat of size R. 

In fig. y we have plotted rigas/Zni at z = 4 computed 




FIG. 1. The neutral gas density Qgas/fui at z = 4 plot- 
ted as a function of a for a CHDM model with f2^ — 0.2. 
The three curved shaded regions are bounded by the Gaus- 
sian PDF (solid line) and a non-Gaussian PDF with A = 0.5 
(dotted line) and correspond to the region in which the ac- 
tual value for the mixed scenario actually lies. The circular 
velocities of the three regions are given by Vc = 75km sec ~^ 
(top), 50kmsec~^ (middle) and 25kmsec~^ (bottom). The 
parallel shaded regions correspond to the observational limits 
if /hi = 1.0 (bottom) and /hi = 0.1 (top). Recall that a = 1 
corresponds to pure adiabatic perturbations and a = to 
strings only. Note that the non-Gaussian PDF gives a lower 
value, albeit only slightly, than the Gaussian case for certain 
values of a. This non-intuitive effect is explained in the text. 



assuming Gaussian fluctuations (solid line) as a function 
of a for models with fie = 0.7, flh = 0.1, fi^ = 0.2 and 
h — 0.5 — those parameters which were suggested to give 
the best fit to the observed galaxy clusterings on large 
scales in ref . 0] . Simulations suggest that realistically one 
might expect /hi ~ 0.1 P5 26|, but if one were wanting 
to be conservative in ruling out a scenario, /hi = 1 can 
be used to compute an absolute upper bound on flgas. 
Ha — 1, that is, just adiabatic fluctuations, then one 
sees that for all the values of Vc used, the upper bounds 
are compatible with the observations; the larger values 
of Vc producing more neutral hydrogen than lower ones. 
However, if /hi = 0.1 none of these are compatible with 
the value of f2gas detected and this can be seen to be true 



for a wide range of cosmological parameters [|lj. Even 
more striking from this flgure is that the observed value 
can be achieved very easily by the inclusion of only a 
very small string component a ~ 0.8, and in fact the 
predictions could be made compatible with even smaller 
values of /hi . 

A simple non-Gaussian distribution which has been 
suggested to represent the PDF in cosmic string mod- 
els pg] is the log-normal distribution given by 



PA{y) 



C 



V2^^ 



exp 



1 



x'^iy) - ^x{y) 



(2) 



log{Cy + B), B = exp(ylV2) and C = y/B^ - B^. This 
has a single parameter A{> 0), with the limiting case 
A — corresponding to a Gaussian distribution. In this 
case, one can deduce that the neutral gas fraction is given 

by 



^gas = /ni^bcrfc 



— ^ log C — 

AV2 \ <JR 



B 



rfc 



A 



2\/2 
(3) 



It was suggested in ref. |2|] that A « 0.17 on cluster 
scales [R ~ 10/i~^Mpc), but on the smaller scales un- 
der consideration here the value is likely to be somewhat 
larger |30|, say A « 0.5. We have also included compu- 
tations of Jlgas for a model using this PDF and A = 0.5 
in fig. n^, to illustrate the effects of non-Gaussianity. In- 
tuitively, one might have expected the effects of intro- 
ducing non-Gaussianity to increase figas, but it appears 
this is not always the case. Using a non-Gaussian PDF 
does enhance the creation of more rare objects, that is 
aii » 5c- However, if aR^ 5c^ as can be the case when 
a is small , then the overall normalization, which allows 
the integration over the whole PDF account for all the 
matter in the universe, is reduced. 




for y > —C/B and zero otherwise, where Ax{y) 



FIG. 2. The same quantities as plotted in fig. hi but for 
fli, = 0.3. Note that none of the models are even compatible 
with observations of figas for a = 1, and that the inclusion 
of the effects of non-Gaussianity is now to increase figas in 
contrast to the case of Sli/ = 0.2. 

Fig. y contains the equivalent information for the larger 
value of rii, = 0.3 (and now fie = 0.6). Clearly, such val- 
ues of fijy are incompatible with the observations if a = 
since even the upper bound on figas is woefully short of 
the observed gas fraction. However, for string-inflation 
admixtures of around 50% each the upper bound on the 
gas fraction becomes within range of the observations if 
one assumes Gaussian fluctuations, albeit with /hi > 0.1. 
The effects of non-Gaussianity are now more intuitive, 
since tr/j << 5c and the DLYA systems can be thought 
of as being rare objects once again. 



Therefore, we have shown that CHDM model can be 
made consistent with high redshift observations, if we are 
prepared to resort to well motivated strings plus inflation 
scenarios. The physical origin of this result is the extra 
fluctuations on small scales, due to the strings, but also 
non-Gaussianity can play an important role. If the DLYA 
systems can be thought of as rare objects {Sc » (Jr), 
then this can enhance their production relative to Gaus- 
sian. But if they are formed as la fluctuations then their 
production is slightly suppressed in non-Gaussian mod- 
els. It is clear that the effects of non-Gaussianity will 
be more prevalent at early times when an is small and 
that less objects will formed later as an increases. This 
has interesting implications for the epoch of galaxy and 
cluster assembly Q. 

To conclude we have suggested that the CHDM sce- 
nario, and more importantly Qm = 1, can be resuscitated 
by the inclusion of a string component since this allows 
structure to be formed earlier than in the pure adiabatic 
case. Apart from the flaw of the CHDM models which we 
have discussed here, there are other observations which 
favour the ACDM model and OCDM models mainly be- 
cause they have fl^ < 1. It would be interesting to 
understand how these observations would be modified 
for the non-Gaussian theories discussed here. It is clear 
that dynamical measures, such as the evolution of X-ray 
cluster abundance, could all be made compatible with 
fini = 1 is a sufficiently non-Gaussian model. Work on 
the observational status of this class of models continues. 
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